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Abstract. We consider a single harmonic oscillator coupled to a bath at zero temperature. As is well 
known, the oscillator then has a higher average energy than that given by its ground state. Here we show 
analytically that for a damping model with arbitrarily discrete distribution of bath modes and damping 
models with continuous distributions of bath modes with cut-off frequencies, this excess energy is less than 
the work needed to couple the system to the bath, therefore, the quantum second law is not violated. On 
the other hand, the second law may be violated for bath modes without cut-off frequencies, which are, 
however, physically unrealistic models. 

PACS. 03.65.Ud Entanglement and quantum nonlocality - 05. 40. -a Fluctuation phenomena, random pro- 
cesses, noise, and Brownian motion - 05. 70. -a Thermodynamics 



1 Introduction 

Thermodynamics originally developed as a purely phe- 
nomenological description of the effects caused by changes 

■ in temperature, pressure, and volume on physical systems 
at the macroscopic scale. At the heart of thermodynam- 
ics there are four well-known laws [T]; the zeroth law al- 
lows us to define temperature scales and thermometers 
while the first law is nothing else than a generalized ex- 

■ pression of the law of energy conservation. The second 
law introduces the concept of thermodynamic entropy, 
which never decreases for an isolated system. The third 
law states that as a system approaches the zero tempera- 
ture, the entropy of the system approaches zero. Later on, 
Boltzmann and his followers created and developed sta- 
tistical thermodynamics by reducing the phenomenolog- 
ically described thermodynamics entirely to the scheme 
of classical statistical mechanics. When quantum mechan- 
ics appeared, the statistical thermodynamics had to take 
into account additional factors offered by quantum me- 
chanics, but the overall structure of thermodynamics, its 
fundamental laws, and its meaning fit for macroscopic sys- 
tems remained unchanged since quantum mechanics was 
believed to play no roles at the macroscopic scale. 

A big challenge for thermodynamics arose with the 
miniaturization of a system under consideration [2]; in 
contrast to common quantum statistical mechanics which 
is intrinsically based on a vanishingly small coupling be- 
tween system and bath, the finite coupling strength be- 
tween them causes some subtleties that must be recog- 
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nized. Recent advances in technology have enabled us to 
experimentally study mesoscopic systems and test various 
fundamental concepts. The field of nano electro- mechanical 
systems (NEMS) especially has emerged with a great po- 
tential, e.g., in quantum limit detection and amplification 
[3,4|, and welcher-Weg ('which-path') interferometry [5]. 
Here, the effects of dissipative environments that are neg- 
ligible in macroscopic resonators become detrimental, and 
the noise is, therefore, a major limiting factor in control 
of NEMS resonators. Theoretically, NEMS resonators can 
be modeled as the simplest form in the scheme of quan- 
tum Brownian motion (see [6] for fundamental aspects of 
quantum Brownian motion). Such a development in var- 
ious fields related to the quantum statistical and meso- 
scopic physics has led to considerable interest in the area 
of quantum and mesoscopic thermodynamics, especially 
with the question raised on the validity of the thermo- 
dynamic laws. Discussions about what is the meaning of 
quantum thermodynamics [2j[7] have started and contin- 
ued up to now. 

The validity of the second law was questioned in the 
scheme of quantum Brownian motion [2] , motivated from 
the observation of the fact that a single harmonic oscil- 
lator coupled to a bath at zero temperature has indeed a 
higher average energy value than the uncoupled harmonic 
oscillator ground state (see also [8]), which could not be 
in accordance with the second law in its Kelvin-Planck 
form [T] that a system operating in contact with a ther- 
mal reservoir cannot produce positive work in its surround- 
ings (cf. for a discussion on the validity of the quantum 
third law in the low temperature limit, see, e.g., Refs. [9], 
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[TP]). However, this argument has been shown to be wrong 
by Ford and O'Connell [UJ; by means of the generalized 
Langevin equation they showed, for the well-known Drude 
model for the spectral density of bath modes, that the ap- 
parent excess energy in the coupled harmonic oscillator, 
however, cannot be used to extract useful work since the 
minimum value of the work to couple the free oscillator to 
a bath takes above and beyond this excess energy, there- 
fore, the second law of thermodynamics is inviolate even 
in the quantum regime (i.e., for cases with non-negligible 
coupling strengths at temperature T = without thermal 
fluctuation). Unfortunately they were unable to explic- 
itly connect their result with its model-independent, deep 
quantum origin, thus the validity of the quantum second 
law for a more general form of the spectral density of bath 
modes, J(u>) would still remain an open question; actu- 
ally, in the experimental study of mesoscopic systems one 
might be able to manipulate the spectral density J(ui), 
to some extent, in his own way. In this paper, we would 
like to discuss the second law for various damping mod- 
els. We will first show the validity of the second law for a 
discrete distribution of bath modes by exactly proving the 
second-law inequality in a simple form obtained from the 
general treatment of the susceptibility (see Sec. \3$). Sub- 
sequently, the inequality will be appropriately applied for 
various continuous distributions of bath modes. It is then 
found that for damping models with cut-off frequencies, 
the second law holds, whereas interestingly, we may have 
its violation for damping models with cut-off frequency- 
free J(w), which are, however, physically unrealistic (see 
Sec. HI . Let us begin with a brief review on the basics of 
the quantum Brownian motion. We will below adopt the 
notations used in \12\. 



where we used p = Mq, and the damping kernel and the 
noise operator are respectively given by 

1 N c 2 

= Ti E C08 ^ f ) ; f (*) = - M ^W + 

1 VI llLi UJ ■ 

3=1 3 1 

]T Cj Ixj (0) cosK t) + ?A sm( Uj t) } . (5) 
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Introducing the spectral density of bath modes as a char- 
acteristic of the bath, 

N c 2 

J H = 7r E^-^-^') J (6) 

3=1 3 3 

we can express the damping kernel as 
2 f°° duo J{lo) 



2 r dw 

7<1) -HJ.T 



cos(w t) 



(7) 



Let us apply the Laplace transform to eq. (@|) with the 
aid of [L31T14] 



£{cos(w,i)}(s) = 



£{sm(ojjt)}{s) = 3 2 
s +w j 



(8) 
(9) 



With s = —iuo + + = — i(w + i + ) we then easily obtain 
q u := £{q(t)}{-iuj + 0+) 

= XH [&-iwM{l + 7 (w)}g(0)+M$(0)l , (10) 



2 Basics and its general treatment 

The quantum Brownian motion in consideration is de- 
scribed by the model Hamiltonian 



H — H s + Hi, + Hgb , 



(1) 



where 

H„ = 



2M 2 QH ' \y2m 3 2 3 3 f 



N 



N _2 
~2 3 



H sb = -qy cj Xj + q 

3=1 3=1 3 



2m,- uj 2 



(2) 



Here, from the hermiticity of Hamiltonian, the coupling 
constants cj are obviously real-valued. Without any loss 
of generality, we assume that 



where the Laplace-transformed damping kernel, the dy- 
namic susceptibility, and the Laplace-transformed noise 
operator are, respectively, given by 



• N r 2 1 

^ = 77 E 



M t-r 1 rrij uj 2 uj 2 — uj 2 ' 

j j j 



M to 2 — uj 2 — iui'yioj) ' 



(11) 
(12) 



N 
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OJ 2 — UJ 2 
3=1 3 



Substituting IUJ into CE]), we get 



3=1 



(13) 



< U>2 < ' ' ' < WJV-1 < tON 



(3) 



By means of the Heisenberg equation of motion for p we 
can derive the quantum Langevin equation 

M 'q + M J* ds 7 (t - s) q(s) + M uj 2 q = |(t) , (4) 



where 



N 2 N 2 N 

j=0 3 = 1 ■> J j'=l 



(14) 
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It is known [15] that the susceptibility x(ui) m lO has 
poles at the normal-mode frequencies of the total system 
H in |U, ±uj k with k = 0, 1, 2, • ■ ■ , N, so that 



*2 



iwfc7(a;fe) 



0. 



(15) 



Here, we might be able to say that a specific k = k 
would represent the "system harmonic oscillator" with the 
normal-mode frequency Qk , uncoupled to the "bath" con- 
sisting of the remaining oscillators with u>k, where k ^ k . 
From eqs. fT2|) . ljT5)l . and l|15p . we have a compact expres- 
sion of the susceptibility, 



N 



1 3=1 



M N 



(16) 



fc=0 
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Without any loss of generality, we here assume that 

wo < iDi < ■ ■ ■ < Copf^i < ujn ■ (17) 

The damping function 7(0;) in the frequency domain 
has, besides eq. ifTTj) . another expression which is suitable 
for the case of a continuous distribution of bath modes; 
from eqs. (7|) and J8]) we obtain 



7M 
7M| 



M JO ir uj' 



M 



f 

Jo 



00 do/ J(n/) 



P 



LJ -\-U) 



(18) 



(19) 



We here used the well-known formula l/(x+i + ) = P(l/x) 
iirS(x) for x = uj' — uj. For the simple Ohmic case Jq{oS) = 
Mj uj with an cj-independent constant 7 , we easily have 
7o(i) = 2 r y a 5{t), and 70 (w) = 7o with a vanishing princi- 
pal (or imaginary) part in lfT9|) . while for the Drude model 
where J<j(w) = M ^ ujuj 2 d j(uj 2 + uj 2 ) with a cut-off fre- 
quency u>d, we have ~fd(t) = lo^d e~ Udt , and 



ld(oj) 



loUJ 2 d 



+ I 



loOJdOJ 



uJd - iu 



(20) 



For a later purpose, it is interesting to compare Dy,(uj) 
in eq. l[T4|) with the denominator of the right hand side in 
(fl6j) . Then, we can easily find that 



JV 
fe=0 



A' 

£< 

3=0 



A' 



7(0); ] 



A 

n< 

3=0 



(21) 



Here, 7(0) = j(t)\ t=0 > in eq. (UJ. From this com- 
parison of the denominators at u> = u>n, we also ob- 
tain D^(u>n) < and so ujn < Similarly, we can 
acquire both D^(uji) < for N odd and D%(uji) > 
for N even, which lead to the fact that ujq < Wi for 
any given AT. Further, we can obtain the relationship, 



D<^(ujj) ■ Dy.(ujj-i) < for any j. Therefore, it is found 
that 

<^o < w i < < • • • < Wjv-i < o)jv_i < uj n < uj n . 

(22) 

By using D^(ujq), we can also show that u>o < ujq < ujn 
(see also [IB]). Within this general treatment of the sus- 
ceptibility, we would like to consider the quantum second 
law below. 



3 General validity of the quantum second 
law (discrete bath modes) 

The energy of the system oscillator H s at zero temperature 
can be calculated by means of the partition function Z = 
Tre-? 6 with = l/k B T as 



Tr H, e 



T=0 



E.(0). 



(23) 



/3->oo 



It is well-known ff\\T7\ that the system-bath entanglement 
induced by the coupling term H s b in ^ leads to the fact 
that the system oscillator H s , initially in a pure state 
(here, its ground state with the minimum energy E g = 
huo/2), is not in the pure state any longer but in a mixed 
state with a fluctuation in energy, and so we actually have 
E s (0) > Eg. It was even discussed in [18] that the energy 
fluctuation measurements can provide entanglement infor- 
mation (cf. for a thermodynamical approach to quantify- 
ing entanglement in bipartite qubit states, see [19]). From 
the fluctuation-dissipation theorem [12], we can also easily 
obtain 



uj 2 ) lmx(uj + i0 + ) . (24) 



MO) = -5- / <M<"o 

27T Jr. 



The factor Imx(w + i + ) can be evaluated from eq. (fl6 
with uj — ► uj + i + . By means of the technique used, e.g. 
in [20], eq. ([24| can be rewritten as 



£*(0) = 



h 1 

2 2?ri 



duj ■ 



G(uj) 



(25) 



where G(uj) = — l/{Mx(w)}, and the integration path 
is a loop around the positive real axis in the complex em- 
plane, consisting of the two branches, (00 + ie, it) and 
(— ie, 00 — ie). Therefore, E S (Q) can be exactly obtained 
in closed form from the residues evaluated at all zeroes of 
G(u>) on the positive real axis. It is also interesting to note 
that the entanglement between any pair of the bath oscil- 
lators Hj = p 2 /2 rrij + m,j uj 2 x] /2 with j = 1, 2, 3, • • • , N 
is induced by the system-bath entanglement and the well- 
known entanglement swapping. £21] As a result, we must 
obtain an excess energy for any j, i.e., (Hj)T=o > hujj/2. 
However, the energy of the total system, (H)t=o = Sfclo ^^fe/2 



is clearly not equivalent to (H s }t=o 
(H s + H h ) T =o. 



(H 



j)T=0 
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The minimum work required to couple a harmonic os- 
cillator at temperature T to a bath at the same temper- 
ature is equivalent to the Helmholtz free energy of the 
coupled total system minus the free energy of the un- 
coupled bath [Tjlllj. The Helmholtz free energy can be 
obtained from the canonical partition function Z s ((3) = 
Tr e-^/Trfc e~ 0kb as F(T) = -k B T\nZ s , where Tr h de- 
notes the partial trace for the bath alone (in the absence 
of a coupling between system and bath, this would exactly 
correspond to the partition function of the system only) . 
By means of the normal-mode frequencies uik the partition 
function can be rewritten as 



where 



Z s ((3) = k ~° " fc ~° 1 TT (26) 

3=1 nj=0 

so that we can easily get, for f3 — > oo, 



N 



2 2\ 



N 



(32) 



fc'=0 

N r 2 



P 



mi ujf { (oji + uJk) 2 {oji - w fe ) 2 



Considering each summand Ak from k = N with keeping 
in mind the frequency relationship in (|22| . we see that 
each of the summand is non-negative and so K > indeed! 
Separately from this result for discrete bath modes, we will 
next discuss the second law for continuous bath modes. 
For doing this job, we will consider a continuation of the 
spectral density J (to) from its original form in ([6]). 



JY 



JY 



m = -£%-£< 



, fe=0 



(27) 



With the aid of eq. pi}, it is evidently found that F(0) > 
hujj/2 for any j = 0, 1, 2, ■ • • , N. Further, we have, from 



F(T) 



da; f(u),T)Tm\— In x(w 



i0 H 



(28) 

where f(u>, T) = ksT ln{2 sinh(Scj/2 /cbT)}. Similarly to 
eq. p5j) , we can obtain an integral form of the free energy 
at T = 0, 



4 The second law for continuous bath modes 

For a discussion of the second law for a continuous distri- 
bution of bath modes, we rewrite eq. I|30p as 



K 



h 

-in 



dUJ -G +{ .) 



(33) 

where the subscripts +/— denote the branches (oo + ie, ie) 
and (—ie, oo — ie), respectively, so that 



G+(lo) := G{u)\ 



, ,2 , ,2 i 



M 



M JO TT Ul> 



1 



(34) 



h 1 / , loG'Iuj) 

F(0) = <f> -V^ 

K ' 2 2m J G(uj) 



(29) 



Here, /(u>,0) = ftu;/2. 

Now, we are in a position to exactly formulate the 
quantum second law within this general treatment; from 
eqs. |25|) and ([29)1 with ffl2|l . we easily find an expression 



K:= F(0)-E s (0) = — 



LU 2 f(m) 

~gVT 



(30) 



and, for the validity of the second law, we have to get 
K > for any TV (the number of the bath oscillators) 
and the limit A*" — > oo. Here, K can exactly be evaluated 
from all residues of the integrand on the positive real axis. 
Substituting ([T6| with x(u) = -l/MG(w) and HQ into 
l(30|). we obtain, after a fairly lengthy evaluation of the 
contour integration (see Appendix[A]for details), the exact 
result 



K 



* N 
8M ^ 



(31) 



andG_(w) := G(w)| 



-iO+ 



G* + (u>). Here, we used G(w) 
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iuj~f(uj) with eq. fl9|) for 7+(a>). Therefore, eq. 



33| easily reduces to 



— Im 

2tt 



R' + (u) 



(35) 



where J2+(w) = — ij+(u). 

First, for the Ohmic case, (7+)o(w) = 7 , which is the 
prototype for damping, we easily obtain K = 0. In fact, 
both (E s )q(0) and Fq(0) have the logarithmic divergence, 
however, the same value, namely, 



(£ fl )o(0) = F (0) 



7T JO 



U) (lJ 2 + U3p) 



( W 2 -^) 2 + (7o^) 2 



(36) 



(see also the discussion in the last paragraphs of Sees. 14.11 
and !4.2p . However, the Ohmic model is not so realistic in 
its strict form because the spectral density of bath modes, 
Jo(w) = M^j bj diverges for large frequencies. We there- 
fore introduce a cut-off frequency uj c which leads to J c (u) 
decaying smoothly to zero for large frequencies u > uj c . 
We will first consider the Drude model, where Jd(w) is 
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polynomially decaying for uj > u> c = ujd, and next a damp- 
ing model with J e {uj) being exponentially decaying for 
uj > lo c = uj e . For these damping models, we will be able 
to show that K > 0. Subsequently, we will also consider 
two different damping models without cut-off frequencies 
u> c ; first, the extended Ohmic models where the spectral 
densities J{uj) diverge polynomially faster than Jq(lo), and 
secondly, the extended Drude models with Jd,n{^) diverg- 
ing faster or more slowly than Jq(uj). Interestingly, we will 
observe K < for some of the cut-off frequency-free damp- 
ing models (see Sees. 14.31 and 14. 4j ). 



4.1 Drude model (d) 

We briefly review the second law in the Drude model 
considered in [11] ; it is convenient to adopt, in place of 
(wo, uid, 7o), the parameters (wo,i7,7) through the rela- 
tions 



2 n 

Wn ~ 1 TT ; UJd 



"o JTTi 



n 



7! 



7. 



(37) 



Substituting eq. §2£^ with (J3T|) into eq. (|T2[). we obtain 
the susceptibility 



4.2 Exponentially decaying model (e) 

We now consider a damping model with J e (uj) = M 7 D ui e~^l u 
which, in the limit u> e — > oo, clearly reduces to Jo(w) for 
small frequencies. Substituting this into eq. J7]), we can 
obtain 

7T 1 + (u} e ty 

Applying the Laplace transform [14] to eq. (J40J) with s = 
— i uj + + , it can be found that 



l e 

■ lo 

i — 

7T 



(41) 



,u>/u 



+ e 



Ei 



(see Appendix [C] for the detailed derivation). By using 
this with E[(y) = —Eo(y) = —e~ y /y, we can easily get 
(R' + ) e (uj) and (G + ) e (w), and then introducing a dimen- 
sionless variable A : 



K, 



Uj/uJ, 
2 IT 2 



e , we arrive at the expression 
A (A) 



Im 



d\ 



(i 



/a (A) 



where 

/i(A) 
/a (A) 



A 2 {e A £i(A) - e~ A Ei(A) + i 7r e - A 
A{e A Si(A) 



2 \2 



(42) 



(43) 
A Ei(A)} 



Xd(w) 
1 



LU + iiOd 



M uj 3 + iudto 2 — (luq + 7 u>d) w — i cjq 
1 w + i (/2 + zi + z 2 ) 
M {lo + if2)(uj + izi)(u> + iz 2 ) 



(38) 



(39) 



where z\ = 7/2 + iwi and Z2 = 7/2 — iwi with wf = 
w 2 , - ( 7 /2) 2 . This gives us (G+) d (u) = -l/{Mxd(u)} 
for eq. |35|) . By means of eq. <(39j) , we can even obtain 
the closed expressions for both (E s )d(0) from ([24]) and 
■Fd(O) from lj28|) . We give the detailed derivation of these 
expressions in Appendix which will also be used in Sec. 
14.41 It has been numerically shown in [11] that (E s )d(0) in 

165]) is actually greater than E g = ^^77^, and F d (Q) 

in (|68| is even greater than the (E s )d(0), i.e., Kd > 0. 



For a later purpose, we will also evaluate Kd explicitly for 
various pairs (u)o,w d ) (see Table[2]in Sec. 14.4]) . 

It is noted that in the limit ujd — > 00 (equivalently, 
SI — > 00), we have — > -2—E g (see Appendix [Bj. From 
the comparison between 7d(o>) and 70 (w) (or, equivalently, 
Jd{w) and Jo(w)), this result would be interpreted as Kq — ► 
-^j^Eg. However, it is misleading; 7d(w) behaves only for 
small frequencies, w <C oJd, like in the Ohmic case, which 
corresponds to jd(t) —* 7o(t) only for large times. Ac- 
tually, 7d(t) = 7o^de~ UJdt with w<j — » 00 does not re- 
duce to 70 (t) = 2j S(t) = lim Wd _ >00 ^= 7 ^7e~^ d 4 . 
For the evaluation of if, however, all frequencies, < 
uj < 00, have to be considered. Therefore, we evidently 
get lim^^oo K d > K = . 



+ i w e 7 A e 



(44) 



We numerically evaluate the integration in 11421) for various 
pairs (7 ,w e ) to show that K e > (see Table[TJ. 

From the fact that 7 e (w) behaves like in the Ohmic 
case for small frequencies uj -C ui e , it is also interesting 
to consider the leading behavior of K e for u e — > 00; from 
eq. (|42"ll we can easily get lim a)e _ > . 00 /f e = - ^r 7^ 0, which 
is also different from lim^^oo Kd in Sec. f4.ll This con- 
firms that these limiting values cannot reveal the Ohmic 
counterpart Kq. 



Table 1. K e /E g for various pairs (j a ,ui e ), where E g — 



% (i.e. 



ljq = 1); lim^ 



Ke/Eg = 7o/7T. 





70 = 0.5 


7o = l 


70 = 2 


jo — 5 


0.5 


0.04225 


0.08186 


0.15604 


0.34038 


1 


0.06130 


0.11838 


0.22117 


0.47348 


5 


0.10600 


0.20348 


0.37899 


0.81614 


10 


0.12131 


0.23326 


0.43819 


0.96224 


50 


0.14414 


0.28018 


0.54302 


1.25567 


80 


0.14789 


0.28896 


0.56377 


1.32020 


00 


0.15915 


0.31831 


0.63662 


1.59155 



4.3 Extended Ohmic models (p) 

Let us consider damping models with J p (u>) = M j Q uj (uj/j ) p 
being polynomially divergent with uj. Clearly, the case of 
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is Ohmic. First, we have J\(w) = Muj 2 . By us- 



,iky 

■Pi 



it5(k) + iP(l/fc), we 
which leads to no well- 



ing the relationship J Q dy 
can easily obtain 71 (t) 
defined 71 (w). This (p = 1) is, therefore, physically not 
acceptable. It is not difficult to show that the cases of p 
being odd are not acceptable. 

Next, we consider the case of p = 2. It can be shown 
that 72 (t) =-£ 6"(t) and (*f + ) 2 (u>) = £ - i « «. By 
using this for eq. (j35j) . we can obtain 



h T r 00 , 2w 2 {-w + i5(o)} 

— Im / — ^ : 7. — 

2n J uj — id uj -\- i fj 



(45) 



where a = 7 D + 2 (5(0) and /3 = Wg j . The integral in l(45[) 
diverges logarithmically. This divergence is, obviously, from 
the fact that both (£3)2(0) and £2(0) diverge logarithmi- 
cally, however, differently from the Ohmic case, (£'3)2(0) ^ 
£ 2 (0). In fact, we find that K 2 = -§ {6(0) +70} x 00 < 0, 
which clearly means that the excess energy, (£.5)2(0), is 
greater than the minimum work (or the work in a re- 
versible process), £2(0), required to couple a system to 
a bath. This violation of the second law in the reversible 
process may be understood to emerge from a large amount 
of the energy offer by the bath with J 2 (uj) diverging with 
uj. The infinite value of £T 2 suggests, however, that this 
model would be strictly unrealistic. 



4.4 Extended Drude models (d,n) 

We now consider a more general class of the spectral den- 
sity than J d (ui), which is 



Jd,n(uj) 



J d (uj) = Mj 



<"5 (^ 2 + ^i) 



(46) 



Let us begin with n being odd. First, n = 1. We then have 
Jd,i(oj) = Mj uj d uj 2 /(uj 2 +uj 2 ), which converges to a non- 
zero constant Mj a ujd for large frequencies. Substituting 
this into eq. (O, we can obtain, after some calculation (see 
Appendix \C\ for details), 



7d,i(*) 



JpUJ d 

ir 



[{Ei(u) d t) + E 1 (u d t)} sinh(uj d t) 



{Ei(ui d t) - E^Udt)} cosh(uj d t) } . (47) 
Applying the Laplace transform [14] to this, we can get 



Here, we introduced a dimensionless variable A = uj/^ 
with Ao = wo/7o and A^ = ujd/j - We numerically evaluate 
Kd t i for various pairs (Ao, Ad) to show that £^,1 > (see 
Table EJ)- It is also noted that in the limit Ad (or u) d ) 
— > 00, the spectral density J d< i(u>) with j = Ud reduces 
to J\((jj) in Sec. 14. 31 for small frequencies. As was discussed, 
however, this case (d, 1) is a well-defined damping model 
whereas the model with Ji(u>) is not. 



Table 2. K d ■K/ r y E g from dTTJ) versus £ d ,i n/-y E g from (J49J) 
for various pairs (ujg,ujd) with 7o = 1, where E g — — j -; (d, 0) 
denotes the Drude model. 



(«0, Ui) 


(d,0) 


(d,l) 


(wo, U>d) 


(d,0) 


(d,l) 


(0.5,0.5) 


0.84275 


0.50184 


(1,0.5) 


1.21124 


0.25468 


(0.5,1) 


0.61942 


0.45203 


(1,1) 


0.45318 


0.26521 


(0.5,5) 


1.29483 


0.34454 


(1,5) 


0.61427 


0.17800 


(0.5,10) 


1.52266 


0.24507 


(1,10) 


0.73767 


0.05454 


(5,0.5) 


0.53089 


0.02832 


(10,0.5) 


0.28968 


0.00936 


(5,1) 


0.47010 


0.04106 


(10,1) 


0.27637 


0.01345 


(5,5) 


0.09790 


0.05476 


(10,5) 


0.13428 


0.02346 


(5, 10) 


0.09537 


0.04976 


(10, 10) 


0.04792 


0.02347 



Let n = 3 next. We have J(w)d,3 = M 70 uj d 1 uj 4 / (uj 2 
I. After a straightforward calculation, we will obtain 



7d,3(*) = -7d,l(*) 



270 p 

ITUJd t 2 



(51) 



which indicates that this case is physically not acceptable. 
Similarly, we can also show that all cases for n > 3 being 
odd are not acceptable. 

Now, let n be even. We can then find that 



7d, 2m (i) = (-l) m { ld(t) 



E 



7o {2(J - 1)} («) 



(52) 



where to = 1, 2, • • • , and 7g 2 ^ 1 ^(i) represent 2(j — 1)- 
time derivatives of 70 (t). We begin with a simple case 



(7+) d i(w) = + i 1 In(^) . (48) (m - 1) with J«j, a (w) = M 7o w 3 /(w 2 + w 2 ). This case is 



Using eq. (|34j) with (|48j) . we can arrive at the expression 
in Qg) 

jr **>t 5i(A) 
£d,i = — Im / dA— — -, 
2tt 7 .92 (A) 



(49) 



where 
5i(A) 



A 2 



2 ^{(A 2 -A 2 ) In (A) + A 2 + A 2 } + 

i\ d (\ 2 -\ 2 d)\ , 
.92(A) = (A 2 + A 2 ) {(A 2 -A 2 )(A 2 + A 2 )- 

^ln^A) + i A d A 2 | . (50) 



particularly interesting because J d> 2(uj) diverges for large 
frequencies, however, more slowly than Jq(u>) for the Ohmic 
case, whereas all Jd.im(u) for to > 1 diverge faster than 
Jo(to); Jd.i(u) may be said to be of weak divergence. Due 
to the fact that £d,o > £d,i > seen from Table [2j 
we would like to pose a question if we will here obtain 
K d ,i > K d .2 > = K or K d;1 > > K d . 2 . In fact, we 
have an interesting relation (^ + ) d , 2 (uj) = —^d(uj) + Jo(uj) 
from eq. |52|) . and so (i?^)^^^) = — (R' + )d(uj) for eq. 
(|35l) . Introducing the parameters (wo, 7) defined as the 
relations 



UJd UJq 



Qvtl ; w d +7o := 



12 + 7; w 2 := l?7 + ' 



(53) 
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(note that these differ from the relations in ij37|) ). we can 
easily obtain 



(G+k 2 M 



uj 3 + i (7o + u>d) uj 2 — WqUJ — iajQOJd 

UJ + lUJd 

(u) + iff) (uj + iz\) (uj + 1Z2) 



uj + ittw 2 /(f2j + w§) 



(54) 



where z\ = ^ + iwi and Z2 — 5 - iwi with w 2 = 
Wq — (^) 2 . By using eqs. ([351 and |54"|). we arrive at the 
expression 



K d , 2 -- 
where 



2n (f> 7 + w 2 ) (f2 7 ~ f> 2 - w 2 ) 



C(w ,/2,7) 
(55) 



C(w , J7, 7) = 7 (w 2 , - I2 2 ) -±- arctan ■ 



Wi 7 

(^ 2 +w 2 - fl 7 ) ln(/2 7 + w 2 ) - 
2(i? 2 + w 2 ) In w + 2^7 lni?. (56) 



In case that Wi is complex- valued (w < 7/2), i.e., for 
the overdamped case, this has to be understood in terms 

where 



1 ■ arctan 2 Wl 



In 



7+2 Wl 



of the relation, 

3 wi 7 iwi \ 7 — 2 wi 

w? = (i) 2 -w 2 . 

Interestingly enough, we can here observe Kd^ < 
(see Fig. [TJ, which would allow us to have a violation of 
the second law in the reversible process for this damp- 
ing model. This negativity may be understood from the 
comparison, with the aid of the relation (R' + )d,2(w) = 
— (R'+)d(w), between 



K d ,2 = T^^dlo x 



(57) 



[uj + itOd) (uj + ifi) (uj + iz±) (uj + IZ2) 



with (wq, Q, 7) in eq. (j37|) and &d = O + 7, and 



K d = -^uJdjo x 



ImJ^duj 



(58) 
> 



with (wo,/2, 7) in lj53|) and = w„ J7/(i7 7 + Wq). In 
fact, we can also show, by using eq. lf28|) with 1(54]). that 
*d,2(0) > F d (0). 

Next, we briefly consider the case of (d, 4). We then 
have (7+)d,4(w) = 7d(w) - 7o + (70 ^ 2 - 2 170 <5(0) cj)/u^ 
from eq. (|52| . After a fairly lengthy calculation with this, 
we can eventually obtain an explicit expression for Kd,4 = 
&f™duf d ,t(w), where f dA (oj) -2 {6(0) + ujj/j^/uj 
for large w. From this asymptotic form, we easily see that 
Kd,A — * —00, which indicates the violation of the second 
law. This infinity of KdA suggests, however, that this case 
would be strictly unrealistic. 

From the above results for cases (d, n) including the 
case (p = 2) in Sec. I4.3[ we may be able to say that aside 
for physically unacceptable damping models, the diver- 
gence (weak or strict) of the spectral density J(uj) for large 
frequencies could lead to the violation of the second law. 



It is also interesting here to note that Re (7+)d,2(w) = 
1o uj 2 /(uj 2 +uj 2 ) > and Re (j+) dA (uj) = j a uj^ 2 uj 4 /(uj 2 + 
u>d) > 0. It is known [23] that a violation of the posi- 
tivity for Re7+(w) is tantamount to a violation of the 
second law in the thermodynamic limit (where a coupling 
strength between system and bath vanishes). We see here, 
however, that this positivity would not be a sufficient con- 
dition for the second law in the quantum regime (with a 
non-negligible finite coupling strength between system and 
bath). 



5 Conclusions 

In summary, we have extensively studied the second law 
in the scheme of quantum Brownian motion. It has been 
observed that from the system-bath entanglement, a sys- 
tem oscillator coupled to a bath at zero temperature has 
a higher average energy value than the ground state of 
an uncoupled harmonic oscillator. For a damping model 
with arbitrarily discrete bath modes and damping mod- 
els with continuous bath modes with cut-off frequencies, 
however, this apparent excess energy has actually been 
found to be less than the minimum work to couple a sys- 
tem to a bath. Therefore, the second law holds in the 
quantum regime. We also found, on the other hand, that 
the violation of the second law may happen for some cut- 
off frequency-free damping models, which are, however, 
physically unrealistic; especially the case (d, 2) in Sec. 14.41 
with a less diverging spectral density of bath modes than 
the Ohmic model being the prototype for damping, has a 
finitely negative value of Kd.i- The further question about 
the validity of the quantum second law for a broader class 
of quantum systems than the quantum Brownian motion 
considered here, particularly non-linear systems coupled 
to a bath, clearly remains open. 
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A : A detailed derivation of eq. ( f3TT ) 

By substituting eq. lfT6|) with = —l/MG(u>) and eq. (fTTj) into K in l(30|). we immediately have 



■t /■ 11 3 N 2 

A- = jp™ I fa J 2 ^1 y —2l — J * + 1 L , (59) 



1 W - u k) 



k=0 

where the integration path over a; is a loop around the positive real axis in the complex w-plane, consisting of the two 
branches, (oo + ie, ie) and (— ie, oo — ie). By using the residues at all poles u = uik of the integrand, we can evaluate 
the contour integration. In doing so, we do not have any residue at lo = u>k' when uji = u>k' ■ Accordingly, we can finally 
obtain 

N 



K = m^ Ak ; Ak = A ™ -a™' (60) 



where 



N 



7=1 

JV 

n &i-*i>) 



N 



A W = Qk ; A k 2) = V Pi " 2 + . 1 . (61) 

fe ^2m, w 2 \ (w, + uj k ) 2 (w, - ^ fc ) 2 ' 



(2) 

Here, ^ is, obviously, positive-valued. Therefore, the non-negativeness of Ak can be completely determined by the 

factor A k . Keeping in mind the frequency relationship in (j22| . we first consider A^ ■ This is clearly non-negative. 
Next, for k = N — 1, we have 

N-i 

^-r - S x "f-i . (62) 



[ (un-i - ^k) 



to 



N-l N 



k=0 



The first factor on the right hand side is non-negative, and so is the second factor whose numerator and denominator 
are negative-valued, respectively. Therefore, we get An-i > 0. Along the same line, we can straightforwardly show 
that each summand Ak with k = N — 2, JV — 3, • • ■ , is non-negative, which will yield K > 0. 



B : No violation of the second law in the Drude model 



In the Drude model, we can even evaluate the system energy (E s )d(0) and the free energy Fd(0) explicitly and in closed 
form. From eq. (|39|) we see that for wo > 7/2 (underdamped case), z\ and Z2 are conjugate complex numbers to each 
other, while for Wo < 7/2 (overdamped case), both z\ and Z2 are real-valued. Therefore, Im {{R' + ) d(us) / (G +) d(u))} in 
1(35]). and thus the explicit expressions of Kd, for both cases would differ from each other in parameters (wo, fl, 7). By 
using eq. (|39f we can obtain 

1 (wg + H 2 -7 2 /2) arccos(7/2w ) - 7W1 ln(/2/w ) 



duj Tmxd(w) = < 



M wi (wg - i7 7 + I? 2 ) 

1 ( 7 2 /2 - wg - Q 2 )/2 ■ In (jgjg) - 7 wi ln(J?/w ) 
M wi (wg - ^7 + J? 2 ) 



for wq > 7/2 



for w < 7/2 



(63) 



/•OO 

/ w 2 ImXd(^) 
Jo 



' 1 (w^ + wgi? 2 - l? 2 7 2 /2) arccos(7/2w ) + l? 2 7Wi ln(/2/w ) 
M wi (wg - J? 7 + I? 2 ) 

1 (l2 2 7 2 - 2w4-2w 2 r2 2 )/4.1n(^fi) +l2 2 7Wl ln(fl/w„) 
M wi (wg - Q 1 + Q 2 ) 



for Wq > 7/2 



(64) 



for Wq < 7/2 
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where wi = \/\Vq — (7/2) 2 and Wi = -iwi. 

From eqs. (|24|) . <[63j) , and ([64|. we have an exact expression 



(E s ) d (0) = — {A(w ,/?, 7 ) + S(w ,r2, 7 )} 



(65) 



where 



' (w 2 + J? 2 )(2J2w 2 + w 2 7 ) - f2 2 7 3 /2 



A(w ,fl,7) 



wi (/2 + 7) (w 
(w2 + fl 2 )(fif/4 



i?7 + I2 2 



arccos(7/2wo) 



w 2 7 /2) - l2 7 2 /2-(w 2 - fl 7 /2) ln ^ 7 /2 



(66) 



wi(/2 + 7 )(wg - ^7 + 
for Wo > 7/2 and Wo < 7/2, respectively, and 

i? 7 (l2 2 + J? 7 - wg) 



Wi 



7/2 



Wi 



B(w ,fl,7) 



(fl + 7)(wg - tf 7 + ^ 2 ) 
Similarly, the free energy fd(0) in eq. (|28|) can be exactly evaluated as 



ln(fl/w ) 



(67) 



-^ ( r2 + 7 ) In 



fl + 7 
Q 

n 



7 In ( — 
w 

7 In ( — 
w 



2wi arccos 
'7/2 



7 



wi In 



2w 

- Wi 



7/2 



Wi 



for w > 7/2 



for wp < 7/2 . 



(68) 



Clearly, (E s )d(0) and Fd(0) for the underdamped case are identical to eqs. (10) and (14) in [TT], respectively. These 
expressions can also be applied for the overdamped case (wi ^ R) with the aid of the complex-valued expression, 
arccos(y) = f + i ln(i y + \J\ — y 2 ) and actually equivalent to those for the overdamped case derived here. Then, we 
can obtain E s (0) < F(0) for both cases (see Fig. 2 in [TT]). 

In the limit fl — ► 00 (equivalently, u>d — > 00), we get from eqs. (f6o|) - <|68|l 



which reduces to 



a; 



7 



ft/2 
"2V 



In 



E„ 



ft •) 



^ + 7 
i? 

1 



2^ 



2QJ 



O 



1 

T? 2 



1 



in) 



o 



1 

J2 2 



7TW " 27T ^1 + -f/n 

(cf. eq. (15) in [HJ). Lastly, we explicitly give an explicit expression for Kd in original parameters (uiq, uj d , 7o ), 

^7o 



A rf = 



2- 



r/A 



2 A 2 A 5 



(2 Aq + Ad) A 2 , A 3 



{(A 2 + A 2 )(A 2 -A 2 )-A d A 2 } 2 + (A 2 A)^ 



(69) 
(70) 

(71) 



where a dimensionless parameter A = ll>/j with Ao = wo/7o and A^ = u) d /~f . Here, we see that the non-zero value of 
K d for Xd — > 00 arises from the competition between two terms of the numerator of the integrand for large A. Eq. (|7T|) 
is also used in Table [2] for comparison with Kd t i in eq. I|49p . 



C : Details for eqs. (SID and (TO 

In derivation of eq. fiTj) . we used the relationship [24| E\(—y ± i0 + ) = — Ei(y) ^ * 7r, where the exponential integrals 
Ai(y) = Jj 00 dz e~ yz /z and Ei(y) = P dz e z /z. Also, for eq. (07]), we employed, first, [25] 

/■°° cosfaw) 

_/ ~ + I = -sin(ab) si(ab) - cos(ab) Ci(ab), (72) 

where the sine integral si(y) = — dz SH^fl = _f _(_ Si(y) with Si(y) = dz sm ^ z ) ) and the cosine integral Ci(y) = 
— J 00 dz £££M = c e +ln y+ dz cos ^)~ 1 with the Euler constant c e = 0.5772156649 • • • ; secondly, we used the relations 
Si(i») = i{Ei(y) + £i(y)} and Cifey) = |{Ei(j/) - ^(y)} + ft [24]. 
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